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CHAPTER  1: 

Introduction 


The  family  of  Ill-Nitride  (III-N)  semiconductor  materials  is  used  extensively  in  electronic 
devices  requiring  high-power  density,  high-temperature  and  high-frequency  operation.  Due 
to  the  wide  band-gap  and  high  mobility  of  III-N  materials,  devices  fabricated  using  III-N 
materials  outperform  devices  fabricated  from  related  III-V  semiconductors  in  breakdown 
voltage,  radio-frequency  (RF)  power  gain  and  thermal  performance.  The  ability  to  fab¬ 
ricate  III-N  heterojunctions  which  lead  to  the  formation  of  a  high  mobility,  high  density 
two-dimensional  electron  gas  (2DEG)  allow  for  the  fabrication  of  high  power  RF  transistors 
for  use  in  applications  where  vacuum  tubes  have  traditionally  been  the  only  option.  Ad¬ 
ditionally,  due  to  the  high  atomic  threshold  displacement  energy  of  III-Ns,  these  materials 
promise  superior  reliability  in  environments  where  radiation  effects  are  a  concern. 

The  2DEG  formed  in  III-N  he tero structures  can  be  leveraged  to  create  a  three-terminal 
field-effect  transistor,  known  as  a  high-electron  mobility  transistor  (HEMT),  which  is  opti¬ 
mally  suited  for  RF  applications.  III-N  HEMTs,  particularly  those  fabricated  using  (0001) 
oriented  Gallium  Nitride  (GaN)-Aluminum  Gallium  Nitride  (AlGaN)  he  tero  structures,  are 
currently  found  in  many  of  our  next  generation  weapons  and  radar  systems.  High  bandwidth 
shipboard  communications,  phased-array  radar  and  electronic  warfare  modules  are  exam¬ 
ples  of  systems  that  are  already  taking  advantage  of  the  properties  wide  band-gap  III-Ns 
provide  [1]— [3].  Due  to  the  increased  power  density  seen  in  GaN  HEMTs,  they  are  effective 
replacements  for  klystrons,  which  take  significantly  more  space  in  shipboard  applications. 
Another  application  for  this  technology  is  satellites.  Again,  because  of  their  decreased  size, 
modules  leveraging  GaN  HEMTs  save  payload  space  without  sacrificing  capability.  This  is 
crucial  for  payload  launch  cost  in  space  systems.  The  effects  of  radiation  upon  electronics 
in  the  space  environment  has  been  well  documented  over  the  years  and  provides  engineers 
with  different  challenges  as  compared  with  terrestrial  environments.  Knowing  how  radi¬ 
ation  affects  device  characteristics  and  reliability  is  of  prime  importance  to  the  design  of 
devices  that  will  be  used  in  space-based  applications.  Extensive  research  upon  the  effects 
of  radiation  upon  the  III-N  material  family  has  been  carried  out,  but  the  understanding  of 
the  effects  in  devices,  especially  in  GaN  HEMTs,  is  not  complete. 
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1.1  Research  Objective 

The  effects  of  non-ionizing  damage  induced  by  high  dose  massive  particle  radiation  on  GaN- 
on-silicon  (GaN-on-Si)  ffEMTs  are  explored  using  a  physics-based  technology  computer- 
aided  design  (TCAD)  model.  The  model  is  based  upon  GaN-on-Si  ffEMTs  fabricated  at  the 
U.S.  Naval  Research  Laboratory  (NRL)  for  radiation  experimentation.  The  Silvaco  ATLAS 
TCAD  package  is  used  to  first  simulate  the  transport  of  charge  carriers  through  the  initial 
GaN-on-Si  ffEMT  structure  based  upon  measurements  taken  at  the  Naval  Postgraduate 
School  (NPS).  Special  attention  is  paid  to  accurate  modeling  of  the  mobility  of  the  2DEG 
channel,  which  is  degraded  by  the  presence  of  radiation  damage. 

Utilizing  the  stopping  and  range  of  ions  in  matter  (SR1M)  software  package,  we  used  a 
Monte  Carlo  simulator  of  heavy  ion  radiation  damage  within  materials  to  generate  damage 
profiles  for  interstitials  and  vacancies  induced  by  proton  irradiation  of  varying  energies  for 
the  ffEMT  structure.  The  SRfM  program  utilizes  the  Kinchin  and  Pease  (K-P)  model  of 
Frenkel  pair  generation  to  simulate  damage  events  under  heavy  ion  irradiation.  Profiles 
were  created  for  proton  energies  of  2.0,  5.0,  10.0,  20.0  and  40.0-MeV  protons.  These 
profiles  were  then  inserted  into  the  ATLAS  model  developed  for  the  initial  device.  A 
literature  review  was  conducted  to  identify  and  classify  the  energy  levels  of  the  related 
traps  associated  with  vacancies  and  interstitials  created  by  the  radiation  damage  predicted 
by  SR1M.  Again,  care  was  taken  to  select  a  physically  realistic  model  of  the  mobility  term 
which  is  affected  by  the  presence  of  radiation  induced  defects. 

Results  were  compared  with  data  taken  from  prior  work  at  NPS  and  NRL  involving  2.0- 
MeV  proton  irradiation  of  the  NRL  devices  to  the  model  outputs.  After  modeling  the  effects 
of  2.0-MeV  proton  irradiation,  we  expanded  the  model  to  predict  damage  characteristics  at 
the  additional  proton  energy  levels,  with  emphasis  on  predicting  the  expected  change  from 
40.0-MeV  proton  irradiation  for  fluence  levels  above  1015  p+  •  cm-2. 


1.2  Related  Work 

Studies  conducted  on  proton  radiation  in  AlGaN/GaN  ffEMTs  often  have  a  wide  variety 
of  design  variation  between  the  devices.  The  manufacturing  process,  materials  used  and 
dimensions  of  the  device  all  play  a  large  role  in  the  characteristics  of  the  device.  While  it  is 
impossible  to  find  independent  research  with  the  identical  devices  utilized,  ffu  et  al.  have 
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done  similar  studies  on  1 .8-MeV  proton  irradiation  of  AlGaN/GaN-on-sapphire  devices  that 
are  applicable  to  the  study  of  the  effects  of  proton  radiation  damage  at  similar  and  increased 
energy  levels  [4]— [6] .  Modeling  of  radiation  effects  has  also  been  studied.  Hui-Chia  Yu 
wrote  an  extensive  paper  on  concepts  of  modeling  diffusion  of  diluted  species  to  varying 
complexities  [7],  [8]. 

At  NPS,  research  on  the  radiation  hardness  of  AlGaN/GaN  HEMTs  has  been  ongoing 
since  the  early  2000s.  Initial  research  was  carried  out  on  commercial  off-the-shelf  devices. 
Thesis  work  at  NPS  started  with  properties  research  and  failure  analysis  before  moving 
to  irradiation  effects  around  six  years  ago  [9].  More  recently,  research  has  been  focused 
on  the  physical  inspection  of  radiation  damage  in  NRL  fabricated  GaN-on-Si  HEMTs 
and  the  irradiation  of  these  devices.  In  conjunction  with  NRL,  work  was  carried  out  to 
irradiate  and  characterize  GaN-on-Si  devices  using  2.0-MeV  protons.  Wade  examined 
the  micro-structure  of  damage  created  due  to  the  irradiation  by  2.0-MeV  protons  at  a 
fluence  of  6.0  x  1014  protons/cm2  [10].  In  that  study,  upon  physical  inspection  of  imagery 
provided  by  scanning  transmission  electron  microscopy,  voiding  was  noticed  under  the  gate 
fingers  of  the  device.  This  voiding  was  not  originally  present  in  the  device.  Along  with 
this  observation,  energy  dispersive  spectroscopy  analysis  preformed  by  Wessel  revealed 
diffusion  of  Nickel  into  the  Gold  layer  [11].  This  diffusion  at  layer  boundaries  was  also 
evident  at  the  AlGaN/GaN  interface.  Fast  neutron  irradiation  experimentation  was  also 
carried  out.  In  an  experiment  by  Iobst,  NRL  GaN-on-Si  HEMTs  were  irradiated  with 
gamma  radiation  to  document  the  effects,  followed  by  irradiation  with  fast  neutrons  (between 
1.0  and  20.0  MeV)  at  the  McClellan  Nuclear  Research  Center.  The  effects  of  the  gamma 
radiation  were  filtered  from  the  effects  of  the  fast  neutron  irradiation  for  analysis  [12]. 

1.3  Thesis  Organization 

This  thesis  is  laid  out  as  follows.  The  device  characterization,  basic  physics  of  proton 
interactions  in  a  crystalline  structure  and  principles  of  vacancy  and  interstitial  creation  and 
diffusion  are  reviewed  in  Chapter  2.  The  integration  of  physical  properties  into  the  Monte 
Carlo  simulation  software  SRIM  along  with  development  of  functions  of  damage  related  to 
depth  based  on  simulation  results  are  discussed  in  Chapter  3.  The  process  for  implementing 
these  functions  in  the  TCAD  software  Silvaco  ATLAS  along  with  the  changes  made  to 
mobility  modeling  in  order  to  more  accurately  encapsulate  all  of  the  effects  contributing 
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to  the  changes  in  mobility  exhibited  in  the  device  are  also  developed.  In  Chapter  4,  the 
implementation  of  the  radiation  effects  into  the  TCAD  model  is  presented.  The  results  of 
the  model  in  comparison  to  physical  results  seen  from  2.0-MeV  proton  radiation  is  analyzed, 
and  predictions  are  made  upon  the  damage  induced  by  higher  energy,  higher  fluence  proton 
irradiation.  Conclusions  and  opportunities  for  further  work  can  be  found  in  Chapter  5. 
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CHAPTER  2: 
Background 


2.1  Properties  of  Gallium  Nitride 

Device  quality  III-N  material  has  been  available  for  device  fabrication  since  the  1990s.  For 
power  and  RF  applications,  the  most  important  material  in  the  III-N  family  is  GaN.  GaN  is  a 
wide  band-gap  semiconductor  with  a  band-gap  of  approximately  3.4  eV  at  300  K,  dependent 
upon  its  crystal  structure.  The  crystal  structure  of  GaN  can  be  either  zincblende  or  wurtzite. 
The  atomic  structure  and  band  diagram  for  wurtzite  GaN  are  shown  in  Figures  2.1  and  2.2, 
respectfully.  Wurtzite’s  crystal  structure  is  hexagonal  and  is  defined  by  two  different  lattice 
constants  a  and  c.  The  wurtzite  crystal  structure  is  different  from  other  crystal  structures 
in  that  it  lacks  inversion  symmetry.  Due  to  this  lack  of  inversion  symmetry,  GaN  (and 
other  III-Ns)  are  pyroelectric:  they  possess  a  net  spontaneous  polarization  Psp  aligned 
antiparallel  along  the  hexagonal  c-axis  of  the  crystal.  Additionally,  III-Ns  are  piezoelectric 
and  modulate  the  polarization  of  the  material  dependent  upon  the  strain  in  the  lattice. 


i  Energy 


Figure  2.1.  Crystal  Structure.  Source:  [13]. 


Figure  2.2.  Band  Diagram.  Source:  [14]. 
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GaN  is  not  an  easy  material  to  grow  in  bulk  at  reasonable  cost,  so  it  is  often  grown  on  a 
substrate.  The  use  of  a  specific  substrate  can  alter  the  material  properties  of  the  GaN.  When 
creating  devices,  GaN  is  not  grown  directly  on  the  substrates  because  differences  in  lattice 
constants  result  in  strains  that  cause  large  defects  in  the  GaN.  This  mismatch  increases  the 
dislocation  density  within  the  active  region  of  the  device,  resulting  in  diminished  device 
performance.  For  that  reason,  buffer  layers  are  grown  on  the  substrate  before  growing  GaN. 
Each  buffer  layer  absorbs  a  portion  of  the  stress  so  that  when  GaN  is  deposited,  the  stress 
is  at  an  optimal  level  for  crystalline  growth.  Traditional  substrates  for  III-N  growth,  such  as 
sapphire  and  silicon  carbide,  can  be  expensive  materials  to  manufacture  and  do  not  have  the 
technological  maturity  of  older  materials  such  as  Si.  The  desirability  of  Si  as  a  substrate 
for  growth  stems  from  the  ability  to  grow  cheap,  large,  high  quality  substrates;  however, 
the  lattice  mismatch  between  Si  and  III-Ns  is  severe.  This  necessitates  the  use  of  several 
additional  buffer  regions  for  successful  III-N  growth.  In  GaN-on-Si  devices  studied  in  prior 
research  at  NPS,  aluminum  nitride  (AIN)  in  combination  with  varying  mole  fractions  of 
AlGaN  were  used  as  transitional  layers  before  the  GaN  layer. 

The  intrinsic  strain  developed  in  a  crystal  lattice  is  the  result  of  an  offset  in  the  lattice 
constant  between  the  growth  layer  and  the  substrate.  This  is  referred  to  as  lattice  mismatch 
strain 

as  ~  aL  , 

e- - ,  (2.1) 

aL 

defined  by  the  lattice  constants  of  the  substrate  as  and  the  lattice  constant  of  the  layer  making 
contact  with  it  cil .  While  too  much  strain  can  be  a  major  issue  during  crystal  growth,  strain 
can  be  engineered  to  be  useful  in  III-N  structures.  Once  a  high-quality  III-N  layer  is  grown, 
additional  layers  of  III-Ns  can  be  grown  pseudomorphically,  in  which  the  entire  epitaxial 
layer  is  under  strain  imposed  via  growth  atop  the  substrate  III-N  layer.  These  additional 
layers  create  high  quality  heterojunctions  which  can  be  used  for  device  applications. 

2.2  Spontaneous  and  Piezoelectric  Polarization 

Polarization  in  III-Ns  manifests  as  both  spontaneous  Psp  and  piezoelectric  Pp-  polarization. 
The  sum  of  these  is  the  total  polarization 


p  -  p  +  p 
1  ~  r  sp  t  r  pZ. 


(2.2) 
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The  units  of  polarization  are  C/cm2.  The  spontaneous  polarization  can  be  calculated 
for  III-N  alloys  by  linearly  interpolating  between  the  known  values  of  the  two  composite 
materials  of  the  alloy.  For  example,  for  AlGaN  of  mole  fraction  jc, 


Psp,AIGaN(x')  —  P sp,AlN%  +  Psp,GaN(  1  •*•)• 


(2.3) 


Piezoelectric  polarization  is  defined  as 

Ppz  =  2  {eu  ~  e-^j  e,  (2.4) 

where  the  value  e;/-  comes  from  the  piezoelectric  coupling  matrix  and  has  units  of  C/cm2. 
The  Cij  terms  come  from  the  stiffness  matrix  of  the  strained  material  and  have  units  of 
GPa  [15]. 

In  real  he tero structures,  dislocations  due  to  imperfect  growth  of  the  III-N  substrate  layer 
relax  the  strain  in  subsequent  layers.  This  leads  to  a  reduction  in  piezoelectric  polarization. 
The  degree  of  relaxation  can  be  quantified  by  the  unitless  coefficient  tj ,  which  varies  between 
0  and  1.  The  relaxation  coefficient  can  be  used  to  calculate  the  reduction  in  polarization  in 
the  he  tero  structure  due  to  strain  relaxation  [15] 


PT  =  Psp  +  (\-ri)Ppz.  (2.5) 

The  net  effect  of  polarization  is  equivalent  to  the  placement  of  a  fixed  sheet  charge  at  the 
endpoints  of  the  material  layer.  In  a  heterojunction,  the  net  fixed  charge  is  determined  by  the 
polarization  difference  between  the  two  materials.  This  fixed  charge  is  depicted  in  Figure 
2.3(b).  The  sign  of  the  charge  at  the  interface  depends  upon  the  growth  direction  of  the 
heterostructure.  In  an  AlGaN/GaN  he  tero  structure  grown  in  the  (0001)  direction,  this  sheet 
charge  is  positive.  In  order  to  maintain  net  charge  neutrality,  electrons  accumulate  to  screen 
out  this  sheet  charge  at  the  interface.  This  is  called  a  2DEG  due  to  the  fact  that  the  large 
electric  field  in  the  GaN  results  in  the  electrons  being  quantum  mechanically  confined  in  a 
triangular  quantum  well.  The  energy  band  diagram  of  such  a  structure  is  shown  in  Figure 
2.3,  where  the  presence  of  the  2DEG  can  be  clearly  seen  where  the  Fermi  level  goes  above 
the  conduction  band  in  the  GaN  layer. 
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Conduction  Band 
Fermi  Level 

Valence  Band 
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Figure  2.3.  Electronic  Structure  of  an  AlGaN-GaN  Heterostructure.  Adapted 
from  [15]. 


The  band  diagram  and  the  sheet  charge  distribution  of  a  representative  AlGaN-GaN 
heterostructure  are  shown  in  Figures  2.3(a)  and  2.3(b),  respectively,  where  eN+DD  is  the 
charge  from  surface  defects  at  the  top  of  the  AlGaN  layer,  -Q7I.(AlGaN)  is  the  polarization 
of  the  AlGaN  at  the  top  of  the  layer,  +Q7r(AlGaN)  and  -Q7r(GaN)  are  the  polarizations  at 
the  AlGaN/GaN  interface,  +Q;r(GaN)  is  the  polarization  at  the  bottom  of  the  GaN  region, 
and  +Qscr  is  the  charge  from  defects  at  the  GaN-substrate  interface. 
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2.2.1  DC  Characteristics  of  GaN  HEMTs 

A  HEMT  is  a  three-terminal  device  based  upon  the  modulation  of  the  conductivity  of  a 
2DEG  channel  in  a  heterostructure.  Such  a  device  is  easily  fabricated  through  the  deposition 
of  a  Schottky  barrier  upon  the  he tero structure  and  two  ohmic  contacts  to  contact  the  2DEG. 
In  an  ungated  he  tero  structure  the  2DEG  cannot  be  modulated  directly.  The  Schottky  gate 
contact  is  used  to  control  the  concentration  of  the  2DEG  and  can  deplete  the  2DEG  given 
enough  applied  gate  voltage. 


The  DC  characteristics  of  a  HEMT  gatestack  can  be  analyzed  through  energy  and 
charge  balance  considerations.  From  such  an  analysis,  the  values  for  the  gate  threshold 
voltage  [16] 

-A  Ec  A  P 

Vth  =  <pb - -  -  — -  (2.6) 

q  Leff 

necessary  to  deplete  the  2DEG  channel  can  be  calculated,  where  <pi,  is  the  barrier  height 
of  the  Schottky  contact,  A Ec  is  the  conduction  band  energy  level  difference  between  the 
AlGaN  and  GaN  layers,  A P  is  the  difference  in  polarization  between  the  AlGaN  and  GaN 
layers,  C'^.  is  the  capacitance  due  to  the  thickness  of  the  AlGaN  layer  and  the  2DEG  offset 
due  to  quantum  mechanical  confinement  of  the  2DEG,  and  q  is  the  electron  charge.  The 
equilibrium  2DEG  concentration  can  be  calculated  from  [16] 


Wldeg  ~ 


q 


(2.7) 


The  DC  drain  current  for  GaN  HEMTs 

1  Ki2\ 

K  -  Vn)V'd,  - -f  I  (2.8) 

where  W  is  the  width  of  the  device,  L  is  the  gate  length,  /./  is  the  mobility,  Ceff  is  the 
capacitance  due  to  the  thickness  of  the  AlGaN  layer  and  the  2DEG  offset  due  to  quantum 
mechanical  confinement  of  the  2DEG,  V's,  V'ds  and  V'ls  are  the  voltages  at  the  points  shown 
in  Figure  2.4,  and  Vtb  is  the  threshold  voltage,  which  describes  the  drain  current  dependence 
upon  the  local  voltages  beneath  the  gate  contact  and  does  not  take  into  account  resistances 
that  are  present  in  the  source  and  drain  access  region  of  the  device. 


/D  =  —  ^  Ceff 
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Based  on  the  diagram  of  the  device  shown  in  Figure  2.4,  Equation  2.8  can  be  modified 
based  on  voltage  loop  calculations  to  include  the  effects  of  the  access  resistances. 


Figure  2.4.  Source  and  Drain  Resistance  Diagram 


The  relationship  between  the  applied  drain  and  gate  voltages  and  Vgs  and  the 
gate-edge  referenced  drain  and  gate  voltages  V'ds  and  V'gs  are  given  by 


Vds  =  IRS  +  IRd  +  V'ds 


(2.9) 


and 

Vgs  =  1RS  +  V'gs.  (2.10) 

From  considerations  of  the  2DEG  equilibrium  concentration  and  electron  mobility,  the 
value  of  the  resistances  in  the  two  access  regions  around  the  gate  can  be  derived  as  [16] 


Rd 


Lgd 

W  <7  /dn^ldeg 


(2.11) 


where  Lgd  is  the  distance  between  the  gate  edge  and  the  drain  edge.  These  parameters  are 
important  to  the  device  and  affect  the  characteristic  curves  of  the  device.  The  resistance 
values  in  these  regions  have  a  large  affect  on  the  device  at  higher  drain  currents. 


2.3  Non-Ionizing  Radiation  Damage  Effects 

When  an  ionized  particle  moves  through  a  crystal  lattice,  it  interacts  with  the  atoms  in  the 
lattice.  These  particles  may  directly  collide  with  atomic  nuclei  physically  or  may  "collide" 
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with  atoms  through  charge  interaction.  The  probability  of  a  high  energy  ion  colliding  with 
atoms  in  the  lattice  is  inversely  proportional  to  the  energy  level  of  the  ion  [17].  When  a 
collision  does  occur,  conservation  of  energy  and  momentum  is  maintained  and  there  can  be 
a  number  of  outcomes.  If  the  particle  has  enough  energy  to  overcome  the  energy  required 
to  displace  the  atom,  known  as  the  threshold  displacement  energy  Eth,  it  knocks  the  lattice 
atom  out  of  place.  This  creates  an  ionized  atom  that  moves  through  the  lattice  with  an 
energy  dependent  upon  the  nature  of  the  collision,  and  a  vacancy  is  left  where  it  was  before 
the  collision.  This  atom  is  known  as  a  primary  knock-on  atom  (PKA).  Depending  on  the 
amount  of  energy  transferred  to  the  PKA,  we  know  it  either  comes  to  rest  or  continues  to 
move.  If  it  comes  to  rest,  it  either  becomes  an  interstitial  atom  or  replaces  a  vacancy  in 
what  is  known  as  a  replacement  collision.  If  the  atom  continues  to  move,  it  may  go  on 
to  collide  with  additional  atoms.  The  atom  continues  to  move  and  collide  until  its  energy 
decreases  below  Et\x,  creating  additional  damage.  Interstitial  and  vacancies  remaining  after 
the  damage  event  are  known  as  Frenkel  pairs.  A  simple  depiction  of  this  process  is  seen  in 
Figure  2.5. 
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Figure  2.5.  Depiction  of  Particle  Collision  in  a  Crystal  Lattice 
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The  radiation  damage  events  are  classified  as  two-particle  scattering  events,  involving 
only  the  incident  ion  and  the  impacted  lattice  atom.  In  lower  energy  collisions,  there  may 
be  more  than  two  particles  involved  in  a  collision  event.  As  the  majority  of  non-ionizing 
radiation  damage  of  concern  in  semiconductors  is  induced  by  high  energy  particles,  this 
work  only  considers  the  effects  of  two-particle  interactions. 

2.4  Non-Ionizing  Radiation  Damage  in  GaN 

Exposing  the  semiconductor  materials  to  non-ionizing  radiation  damage  causes  several 
different  effects.  The  first  is  damage  to  the  lattice.  When  atoms  are  displaced  from  their 
ordered  lattice  locations  to  a  more  disordered  location,  electron  trap  states  are  created. 
There  are  several  types  of  trap  states.  They  can  be  either  acceptor  or  donor  traps.  Both 
types  can  potentially  accept  or  donate  more  than  one  electron  or  hole.  If  a  trap  state  is  filled, 
it  is  said  to  be  ionized.  If  enough  trap  states  are  ionized  in  a  region,  this  results  in  charge 
build  up,  which  can  alter  the  field  levels  within  a  device  and  lead  to  shifts  in  the  device 
characteristics.  This  can  lead  to  undesired  changes  in  circuit  operation. 

In  GaN  layers,  three  trap  states  are  known  to  be  common  after  proton  irradiation  of 
fluence  levels  above  1012  protons/cm2.  These  are  referred  to  as  ER1,  ER2  and  ER3.  The 
"E"  stands  for  electron  traps,  and  the  "R"  indicates  that  it  was  developed  from  high-energy 
particle  irradiation.  Work  by  Auret  in  [1 8]  has  found  that  these  three  energy  levels  are  found 
at  0.13, 0.16  and  0.2  eV  below  the  conduction  band  of  GaN  after  2.0-MeV  proton  irradiation. 
Auret  does  not  specifically  attribute  these  defects  to  specific  defect  types  or  atoms,  and  they 
are  all  shallow  traps.  In  other  work  by  von  Bardeleben,  a  nitrogen  split  interstitial  defect  is 
noted  in  GaN  after  12.0-MeV  proton  irradiation  to  fluence  levels  of  8  x  1016  protons/cm2 
which  occurs  around  1.0  eV  below  the  conduction  band  by  using  electron  paramagnetic 
resonance  [19].  Most  notably,  a  review  of  radiation  effects  by  Pearton  [20]  examined  and 
attempted  to  fingerprint  individual  radiation  induced  defects  in  the  GaN  lattice  and  the  trap 
energy  levels  associated  with  them.  In  general,  radiation  induces  nitrogen  vacancies  which 
act  as  shallow  donors  at  0.06  eV  below  the  conduction  band,  nitrogen  interstitials  acting 
as  deep  acceptors  around  1 .0  eV  below  the  conduction  band,  Ga  vacancies  forming  deep 
acceptors  at  1.0  eV  above  the  valence  band,  and  Ga  interstitials  forming  deep  donors  at  0.8 
eV  below  the  conduction  band. 
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The  crystal  lattice  is  constantly  vibrating  and  moving  because  of  thermal  energy 
maintained  by  the  atoms  and  their  bonds.  This  allows  atoms  in  the  lattice  to  change  their 
position  without  any  outside  influence.  If  the  temperature  is  not  to  the  point  of  melting 
the  material,  the  device  can  anneal  out  vacancies  and  interstitials.  After  irradiation  events, 
devices  try  to  return  to  their  original  form,  and  device  characteristic  have  been  seen  to 
recover  from  these  effects  even  at  temperatures  as  low  as  100  K  [21].  In  this  work.  He,  Ar, 
O,  Au,  Te  and  Xe  ions  were  implanted  in  GaN  grown  on  sapphire  at  15  K  and  Rutherford 
Backscattering  analysis  was  performed  to  identify  movement  of  ions  as  temperature  was 
increased. 

Sometimes,  when  radiation  displaces  an  atom,  the  atom  is  transferred  enough  energy  to 
move  it  substantially  from  its  vacancy.  This  makes  it  difficult  for  it  to  return  and  anneal  out 
of  existence  or  return  to  ordered  form.  When  larger  amounts  of  damage  are  created,  piles 
of  interstitials  and  vacancies  start  to  form,  and  irreversible  voiding  occurs  in  the  material. 
This  has  been  seen  in  work  at  NPS  in  [10]  and  work  by  Limpijumnong  in  [22].  In  [10], 
Wade  visually  inspected  a  HEMT  that  had  undergone  2.0-MeV  proton  irradiation  and  found 
an  accumulation  of  vacancies  manifesting  as  a  void  under  the  Ni  gate  finger  of  the  device. 
Limpijumnong ’s  work  sought  to  classify  and  explain  native  defect  diffusion  in  GaN  crystal 
lattice  through  density  functional  theory. 
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CHAPTER  3: 
Model  Development 


3.1  Device  Composition 

The  GaN-on-Si  wafers  for  the  HEMT  used  in  this  study  were  grown  by  Nitronex  Inc.  The 
substrate  is  350-//m  thick  Si.  A  silicon  nitride  ( SiNx )  interface  layer  with  a  width  of  2.0-nm 
and  an  AIN  buffer  layer  of  300.0  nm  were  grown  as  an  initial  growth  buffer  layer,  followed 
by  two  layers  of  AlGaN.  The  first  is  400.0-nm  thick  with  an  A1  mole  fraction  of  0.7, 
while  the  second  is  300.0  nm  thick  with  an  A1  mole  fraction  0.4.  A  0.8  /im  high  quality 
GaN  layer  was  grown  on  top  of  the  compliance  layer.  A  17.5-nm  thick  AlGaN  layer  with 
an  A1  mole  fraction  of  0.27  was  grown  pseudomorphically  as  a  barrier  layer  to  form  the 
2DEG  channel.  HEMT  Device  fabrication  was  carried  out  at  NRL.  The  device  gatestack 
consists  of  a  30.0-nm  thick  Ni  gate,  followed  by  a  gold  overlay  of  300.0-nm.  The  device  is 
passivated  with  a  100.0-nm  thick  SiNx  layer.  Ohmic  contacts  to  the  2DEG  were  formed  via 
the  deposition  of  a  Ti/Al/Ni/Au  alloy  followed  by  rapid  thermal  annealing.  The  diagram  of 
the  HEMT  gate  stack  structure  is  shown  in  Figure  3.1. 


AI0,Ga03N 


Si  (111)  Substrate 


Alo.27^a0.73^ 


Si3N4 


Figure  3.1.  Representation  of  the  GaN  High-Electron  Mobility  Transistor 
Gate  Stack  Structure 
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3.2  Initial  Device  Characterization 

The  initial  device  electrical  characterization  was  carried  out  using  a  probe  station  and  a 
Keysight  B1500A  semiconductor  parameter  analyzer.  DC  IV  measurements  were  per¬ 
formed  to  characterize  the  I(is  versus  Vgs  transfer  characteristics,  Igs  versus  Vgs  gate  leakage, 
and  Ids  versus  V/v  characteristics  of  the  device.  Thirty-six  devices  on  three  identical  reticles 
were  tested  to  identify  a  device  that  could  serve  as  a  baseline  for  the  development  of  the 
non-irradiated  portion  of  the  model.  The  section  of  the  reticles  considered  had  geometries 
identical  to  those  used  in  2.0-MeV  irradiation  experiments  performed  upon  these  devices 
in  past  work.  The  device  dimension  of  the  selected  device  are  illustrated  in  Figure  3.2  in 
units  of  //m. 


Figure  3.2.  Scanning  Transmission  Electron  Microscopy  Photograph  of  a 
Single  High-Electron  Mobility  Transistor  from  Reticle.  Adapted  from  [11]. 


3.2.1  Initial  Silvaco  Model  Development 

The  initial  model  of  the  device  was  created  in  Silvaco  ATLAS  based  on  the  physical 
dimensions  of  the  actual  device  as  seen  in  Figure  3.3. 
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Figure  3.3.  Device  Created  In  ATLAS  versus  Scanning  Transmission  Electron 
Microscopy  Imaging  of  Device.  Adapted  from  [11]. 


The  mesh  was  emphasized  around  the  GaN  and  AIGaN  interface  and  vertically  under 
the  edges  of  gate.  Initially,  a  built-in  polarization  statement,  polarization,  was  used 
for  the  AIGaN  and  GaN  layers  to  simulate  the  spontaneous  and  piezoelectric  polarization. 
Neither  the  AIN  nor  the  AIGaN  buffer  layers  received  polarization  statements  because  of 
the  defective  nature  of  these  layers  in  the  device.  These  layers  are  assumed  not  to  provide 
any  net  polarization.  Using  these  basic  model  parameters,  we  extracted  the  band  diagram 
seen  in  Figure  3.4,  showing  the  presence  of  a  2DEG  at  equilibrium. 
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AIGaN/GaN  Band  Diagram 


Figure  3.4.  Silvaco  Generated  Band  Diagram  Beneath  the  Gate  Stack 


The  polarization  statement  was  replaced  with  a  sheet  charge  statement  in  the  final 
version  of  the  model.  This  provides  the  same  effect  and  does  not  alter  the  models  charac¬ 
teristics.  The  relaxation  coefficient  rt j,  defined  in  Equation  2.5,  was  adjusted  to  calculate  the 
appropriate  sheet  charge  to  fit  2DEG  density  measured  via  Hall  measurements  [23].  In  or¬ 
der  to  model  heating  throughout  the  device,  a  temperature  dependent  thermal  conductivity 
model  was  utilized,  implemented  via  TCON .  POWER  from  ATLAS’  GIGA  self-heating  mod¬ 
ule,  for  each  material  in  the  device.  The  characteristic  equation  for  the  thermal  conductivity 
model  is 


k(T)  = 


a 

(Tl/300)P  ' 


(3.1) 


The  units  for  k  are  W/cmK  and  the  coefficients  a  and  (3  for  their  respective  layers  are  found 
in  Table  3.1. 


18 


Table  3.1.  Thermal  Heating  Coefficients 


Material 

a 

P 

AIq^iGcIqj^N 

1.381 

0.317 

GaN 

0.65 

0.28 

AIq^GciqjN 

0.688 

1.5535 

AIN 

1.425 

1.64 

Si 

1.48 

0 

Au 

2.0 

0 

Ni 

1.58 

0 

The  source,  drain,  and  gate  were  added  to  the  simulation  with  electrode  statements. 
The  source  and  drain  were  modeled  as  ohmic  contacts  with  work  functions  of  3.6  eV 
to  match  the  affinity  of  the  AlGaN  and  with  an  intrinsic  10.0-D  resistance.  The  10.0-D 
resistance  was  added  to  appropriately  simulate  the  measured  contact  resistance  of  the  ohmic 
contacts.  The  gate  was  modeled  as  a  Schottky  contact  with  a  barrier  height  of  1.1  eV.  The 
Shockley  Reed  Hall  model  was  used  to  simulate  carrier  recombination  in  the  model. 

3.2.2  Mobility  Modeling 

Using  appropriate  mobility  models  is  a  key  aspect  in  accurately  modeling  III-N  devices.  An 
initial,  static  value  of  1260  cm2/(V  •  s)  was  adopted  from  Hall  measurements  taken  on  the 
bare  reticle  [23].  We  quickly  noted  that,  based  on  measured  I-V  characteristics,  this  value 
resulted  in  a  large  overestimate  of  device  currents.  Comparison  between  measured  device 
characteristics  and  the  results  of  this  initial  mobility  model  are  shown  in  Figures  3.5  and 
3.6. 
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Figure  3.5.  Vg  versus  1,j  Vj  =  100.0  mV 


Figure  3.6.  Vg  versus  Id  Vd  =  10.0  V 
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Initial  attempts  at  accurate  mobility  modeling  were  undertaken  by  implementing  a 
custom  mobility  model  within  the  GaN  layer  of  the  device.  ATLAS  allows  users  to  create 
a  custom  mobility  model  which  can  be  implemented  via  a  c-function  through  the  tofimun 
flag  in  the  mobility  statement.  As  a  starting  point,  the  well-known  Albrecht  model 
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(3.2) 


was  utilized  for  GaN  low-field  mobility  [24].  In  the  Albrecht  model,  /u(N,Tl)  is  the 
mobility  as  a  function  of  doping  and  lattice  temperature,  N  is  the  total  doping  concentration 
in  cm-3,  No  is  the  initial  doping  concentration  in  cm-3  and  Tl  is  the  lattice  temperature  in 
K.  Parameters  A,  B  and  C  are  scaling  factors  to  adjust  the  impact  of  the  individual  terms 
of  Equation  3.2.  It  is  important  to  note  that  Equation  3.2  was  developed  for  bulk  GaN. 
Because  the  Albrecht  model  was  derived  for  bulk  material,  it  is  technically  incorrect  to  use 
it  to  model  electron  mobility  in  the  2DEG.  The  first  term  of  the  model, 


A  ■  N 


(3.3) 


accurately  models  impurity  scattering  in  the  bulk  but  not  in  the  2DEG.  The  second  and  third 
term,  used  to  model  acoustic  and  polar  optical  phonon  scattering,  can  be  shown  to  have  the 
same  temperature  dependence  as  the  corresponding  mobility  terms  in  the  2DEG-specific 
model  [25];  therefore,  those  terms  were  incorporated  as  is  and  provided  the  foundation  for 
the  model.  Charge  dislocation  and  alloy  scattering  are  additional  factors  that  this  model 
does  not  individually  tailor  terms  for  but  can  be  shown  as  well  to  be  constants  independent 
of  temperature  and,  thus,  are  incorporated  in  the  tuning  of  the  Albrecht  model’s  B  and  C 
terms.  The  values  for  the  B  and  C  coefficients  were  initially  scaled  to  give  an  initial  mobility 
matching  that  of  the  measurements  in  [26] . 


For  high  field  mobility,  the  Barnes  model  was  utilized  [27].  The  Barnes  model  is 


BnO  +  (Vsat/E)(E/Ecrity 
i  +  (E/Ecrity 


(3.4) 
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where  fi„o  is  the  low  field  mobility,  Vsat  is  the  saturation  velocity,  E  is  the  electrical  field,  Ecrit 
is  the  critical  electrical  field  value,  and  y  is  a  unit  less  exponential  fitting  term.  The  Barnes 
model  derives  from  measurements  of  the  2DEG  electron  velocity  field  scattering  [27]. 
Values  used  for  Vsat  and  Ecrit  were  3.0  x  107  and  20.0  x  103,  respectively.  Based  on  recent 
work,  other  high  field  models  showing  velocity  overshoot  under  conditions  of  high  electric 
field  may  not  be  accurate  [28].  In  the  work  by  Bajaj,  high  fields  were  utilized  to  try  to  reach 
velocity  overshoot  but  found  that  electron  velocity  did  not  overshoot  but  merely  saturated. 
More  importantly,  mobility  overshoot  was  found  to  occur.  This  model  also  exhibits  mobility 
overshoot,  which  fits  the  physical  characteristics  more  accurately. 

Although  modulation  of  the  Albrecht  scale  factors  brought  the  simulated  output  cur¬ 
rents  into  line  with  measured  values,  significant  difficulty  was  found  matching  simulated 
to  measured  IV  characteristics  in  regions  of  bias  where  the  current  was  limited  by  access 
resistances  as  seen  in  Figures  3.7  and  3.8.  To  mediate  this  difficulty,  measurements  were 
taken  to  determine  the  access  resistance  for  the  gate  and  source  regions  of  the  device  and 
can  be  seen  in  Figure  3.9. 


Figure  3.7.  Vg  versus  /,/  Vj  =  100.0  mV 
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Drain  Current  (A) 


Figure  3.8.  Vg  versus  Ij  Vj  =  10.0  V 


Gate  Current  (A) 

Figure  3.9.  Source  and  Drain  Resistance  Measurements 
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From  Equation  2.11,  the  width  in  the  source  and  drain  regions  along  with  the  nideg 
values  should  be  the  constant  in  a  device.  This  leaves  the  mobility  term  and  the  length  of 
the  region  that  differs  between  the  source  and  drain.  If  the  ratio  of  the  resistances  is  not 
equal  to  the  ratio  of  the  access  region  lengths,  then  the  mobility  term  must  be  the  source  of 
variation.  The  ratio  of  resistances  in  this  device  is  280/90,  which  is  not  equal  to  the  ratio  of 
the  access  region  lengths,  13/2. 

From  this  information,  an  additional  factor  which  could  explain  the  measured  deviation 
of  access  region  mobility  from  expected  values  was  sought.  Recent  work  by  Lin  to  explain 
the  discrepancy  in  [29],  [30],  proposed  that  the  mobility  in  the  access  regions  is  reduced  due 
to  an  effect  known  as  polarized  coulomb  field  (PCF)  scattering.  PCF  scattering  is  due  to  the 
fact  that  strain  in  the  AlGaN-GaN  interface  region  is  not  constant  throughout  the  device. 
The  strain  is  changed  by  the  deposition  of  materials  such  as  the  source,  gate  and  drain 
metal  and  the  passivation  layer.  This  alters  the  electron  mobility  as  a  function  of  distance 
between  the  source  and  drain  of  the  device.  As  a  result,  the  mobility  model  was  expanded 
to  include  this  effect.  In  order  to  properly  account  for  PCF  scattering  potentially  induced 
by  deposition  of  additional  layers  above  the  AlGaN,  the  Silvaco  model  was  shifted  from 
ATLAS  to  Athena  in  order  to  utilize  Silvaco’s  Victory  Stress  functionality.  After  running 
the  model  in  Victory  Stress,  non-uniform  strain  under  the  nitride  layer  was  apparent  as  seen 
in  Figure  3.10  which  shows  the  XX  component  of  strain  in  the  lattice.  The  variation  of 
the  strain  causes  scattering  and  a  decrease  in  electron  mobility  from  the  measured  values 
predicted  by  Hall  measurements  on  a  bare  heterojunction. 
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Figure  3.10.  Athena  Victory  Stress  Simulation 


Although  Victory  Stress  could  not  be  used  to  directly  predict  the  effects  of  inhomo¬ 
geneous  strain  in  the  AlGaN  layer  upon  the  mobility,  the  effects  of  PCF  scattering  were 
implemented  in  another  way.  In  order  to  take  into  account  the  PCF  scattering,  the  GaN 
region  of  the  device  was  split  into  three  sections.  This  allowed  us  to  construct  slightly 
different  mobility  terms  for  each  region  based  on  the  resistances  derived  for  each  region. 
The  resulting  constants  used  to  fit  the  mobility  are  presented  in  Table  3.2 


Table  3.2.  Constant  Values  for  GaN  Regions 


GaN  Region  0 

GaN  Region  1 

GaN  Region  2 

A 

0.25 

0.25 

0.25 

B 

4.0 

1.35 

2.0 

C 

4.0 

1.35 

2.0 
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3.2.3  Final  Model  Evaluation 

Our  initial  comparison  between  the  physical  device  and  the  simulated  device  did  not  include 
any  radiation  effects  and  was  used  as  a  measure  for  the  accuracy  of  our  starting  model.  After 
scaling  the  model,  a  strong  correlation  to  the  physical  device  is  shown  in  the  characteristic 
IV  curves  illustrated  in  Figures  3.1 1,  3.12  and  3.13. 


Figure  3.1 1.  Vg  versus  Id  Comparison  Between  Model  Output  and  Measured 
Results,  Vd  =  100.0  mV 
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Figure  3.12.  Va  versus  Id  Comparison  Between  Model  Output  and  Measured 
Results,  Vd  -  10.0  V 


Figure  3.13.  Id  versus  Vd  Comparison  Between  Model  Output  and  Measured 
Results,  Vg  =  -3.0,  -2.0,  -1.0  and  0.0  V 
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3.3  SRIM  Model 

In  order  to  model  the  proper  amount  of  radiation  damage  in  the  AlGaN/GaN-on-Si  HEMT, 
the  Monte  Carlo  ion  radiation  damage  simulator  SRIM  was  utilized  .  A  layered  device  that 
has  the  equivalent  dimensions  of  a  two-dimensional  slice  of  the  HEMT  underneath  the  gate 
was  designed.  These  dimensions  are  found  in  Table  3.3. 

Table  3.3.  Device  Dimensions  Listed  from  Top  to  Bottom 


Material 

Thickness 

S13N4 

100.0  nm 

Au  Gate 

300.0  nm 

Ni  Gate 

30.0  nm 

AIq.2iGcIq.-kN 

17.5  nm 

GaN 

800.0  nm 

AIq.4GcIq.qN 

300.0  nm 

AIq.iGciqkN 

400.0  nm 

AIN 

300.0  nm 

5/3A4 

2.0  nm 

57(111)  Substrate 

350.0  /vm 

After  the  appropriate  dimensions  and  material  names  were  placed  in  SRIM,  the  correct 
densities  were  added.  For  layers  that  consisted  of  a  monochromatic  material,  the  default 
density  was  used.  For  layers  that  have  specific  molar  compositions,  the  density  was  calcu¬ 
lated.  The  compound  correlation  value  was  left  at  one.  The  symbols,  name,  atomic  number, 
weight,  atom  stoichiometry,  Damage  Latt  and  Surf  were  all  left  at  their  default  values.  The 
displacement  threshold  energy  values  found  by  Pearton  were  used  [31].  It  is  noted  that  while 
the  displacement  threshold  energy  levels  for  the  same  atom  vary  in  different  compounds, 
the  assumption  is  made  that  in  GaN,  AlGaN,  and  AIN  the  thresholds  were  relatively  similar 
and,  therefore,  were  left  as  the  same.  The  material  properties  that  were  changed  can  be 
found  in  Tables  3.4  and  3.5.  The  type  of  TRIM  calculation  used  was  a  "Detailed  Calculation 
with  Full  Damage  Cascades."  The  ions  used  to  simulate  p+  irradiation  were  hydrogen  atoms 
with  the  default  mass  and  charge,  zero  angle  of  incidence,  and  2.0-MeV  energy  level. 
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Table  3.4.  TRIM  Compound  Setup  Values 


Compound 

Density  (g/cm3) 

S13N4 

3.2 

Au 

19.311 

Ni 

8.8955 

AIoijGciojjN 

3.0327 

GaN 

6.5 

AIoaGcio^N 

2.8246 

AIqjGcIq^N 

2.3443 

AIN 

3.26 

S73V4 

3.44 

Si{  111) 

2.3212 

Densities  of  composite  materials  were  linearly  interpolated  based  on  molar  composi¬ 
tion. 


Table  3.5.  TRIM  Atomic  Setup  Values 


Atoms 

Displacement 

Threshold 

Energy  (eV) 

Si 

15.0 

N 

28.0 

Au 

25.0 

Ni 

25.0 

Al 

72.0 

Gci 

72.0 

N 

32.4 

Si{  111) 

15.0 

Displacement  threshold  values  from  [31] 
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The  SRIM  program  only  collects  data  in  100  bins  at  a  time  regardless  of  the  size  of  the 
device  being  modeled.  The  larger  the  structure,  the  wider  the  bin  is  required  to  be,  reducing 
the  resolution  of  the  data  collected.  Due  to  the  size  variation  of  layers  within  the  device, 
several  simulations  were  run  with  windows  of  varying  widths.  Emphasis  was  placed  on 
boundary  regions  and  areas  where  current  flows  within  the  device.  The  following  window 
sizes  were  utilized  in  our  simulations:  0-3500,  3500-6500,  6500-22475  and  22475-3522500 
angstroms.  The  window  size  has  no  effect  on  the  simulation  calculations,  it  merely  allows 
the  bin  size  to  be  manipulated. 

Additional  simulations  were  run  at  energy  levels  of  5.0  MeV,  10.0  MeV,  20.0  MeV  and 
40.0  MeV.  Due  to  the  decreasing  amount  of  particle  interaction  with  the  lattice  as  energy 
level  increases,  the  number  of  ions  simulated  increased  for  increasing  energy  levels.  In  order 
to  get  accurate  information  about  the  interaction  of  2.0-MeV  proton  radiation,  10  million 
particles  were  simulated.  For  40.0-MeV  protons,  50  million  particles  were  simulated. 

The  resultant  damage  through  irradiation  was  estimated  through  measurements  of 
interstitials  and  vacancies  from  the  materials  in  each  window.  Energy  levels  less  than  20.0 
MeV  showed  a  linear  trend  in  the  interstitial  and  vacancy  formation  in  both  the  AlGaN  and 
GaN  regions.  The  trend  was  directly  correlated  to  the  thickness  of  the  layer.  The  linear  fit 
for  the  AlGaN  and  GaN  regions  at  2.0-MeV  protons  are  shown  in  Tables  3.6,  3.7  and  the 
interstitial  distribution  in  the  AlGaN  and  GaN  regions  in  Figures  3.14  and  3.15.  The  20.0 
and  40.0-MeV  energy  levels  showed  a  constant  trend  of  damage  in  the  lattice.  The  fitting 
information  for  the  additional  energy  levels  is  found  in  Appendix,  Section  A.  1 . 
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Interstitials  per  Ion  per  cm2  (Int/Ion/cm2) 


Figure  3.14.  Linear  Fitting  for  2.0-MeV  Protons  in  GaN. 


Table  3.6.  Interstitial  and  Vacancy  Distribution  Fitting  Functions  in  the  GaN 
Region  of  the  Device 


Lines  for  2.0-MeV  energy  levels. 


Atom 

Slope 

F-intercept 

Gaj 

1460.437 

-408.6899 

Gay 

1439.269 

-403.156 

N, 

1908.567 

-504.769 

Ny 

1883.116 

-500.607 
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Interstitials  per  Ion  per  cm2  (Int/Ion/cm2) 


Figure  3.15.  Linear  Fitting  for  2.0-MeV  Protons  in  AlGaN 


Table  3 .7 .  Interstitial  and  Vacancy  Distribution  Fitting  Functions  in  the  AlGaN 
Region  of  the  Device 


Lines  for  2.0-MeV  energy  levels. 


Atom 

Slope 

F-i  intercept 

All 

5015.562 

-2143.11 

Aly 

4723.538 

-2015.279 

Gai 

18310.142 

-7809.917 

Gay 

18157.719 

-7742.635 

Ni 

31059.048 

-13205.667 

Ny 

30677.914 

-13042.54 
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CHAPTER  4: 
Model  Implementation 


4.1  Impurity  Scattering  Mobility  Implementation 

The  final  piece  of  the  mobility  term  developed  in  this  work  was  for  the  effect  of  impurity 
scattering.  As  stated,  our  device  was  not  intentionally  doped,  so  an  impurity  scattering 
term  that  is  dependent  on  concentration  of  dopants  returns  a  value  that  has  no  effect  on  the 
initial  model.  Because  of  this,  the  tuning  of  the  model  to  match  the  characteristics  seen  in 
the  non-irradiated  device  was  carried  out  without  an  impurity  dependent  term. 


To  incorporate  the  effects  of  impurity  scattering  in  the  absence  of  the  first  term  of 
Albrecht’s  mobility  equation,  the  Brooks-Herring  model, 


1  afiTi 

HBH  N imp 


(4.1) 


was  used,  where  a  and  / 3  are  user  defined  scaling  factors,  Njmp  is  the  impurity  concentration 
in  cm-3  and  T  is  the  temperature  in  K.  The  Brooks  Herring  model  can  be  shown  to  provide  a 
more  accurate  depiction  of  the  impurity  dependence  of  the  scattering  in  a  GaN  2DEG  [25]; 
however,  the  true  term  should  also  depend  upon  the  2DEG  density  «2 deg>  which  is  not 
included  for  reasons  of  convergence.  The  a  and  /3  terms  were  tuned  to  be  equal  to  0.4  and 
1.0  so  the  model  would  match  the  measured  mobility  percent  change  values  in  [23].  The 
result  of  this  fitting  is  shown  in  Figure  4.1. 
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Figure  4.1.  Percent  Change  in  ji2deg  versus  Fluence 


The  fit  for  the  mobility  is  almost  exact,  while  the  simulation  shows  an  over  estimate 
for  all  the  other  properties  this  study  reviewed. 

4.2  Model  Evaluation  Versus  2.0-MeV  Experimental  Re¬ 
sults 

The  data  on  the  change  in  characteristics  of  the  HEMT  irradiated  with  2.0-MeV  protons 
were  performed  in  [26].  In  [26]  sheet  resistance  Rsh,  2DEG  concentration  n2DEG,  electron 
mobility  n 2deg>  on  state  resistance  Ron,  transconductance  GmjSat,  threshold  voltage  V*  were 
measured  versus  fluence.  The  experimental  results  are  shown  in  Figure  4.2. 
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Radiation  Dose  (cm'  ) 


Figure  4.2.  Percent  Change  Values  after  2.0-MeV  Proton  Irradiation.  Source: 
[26], 


We  matched  the  fluence  levels  in  [26]  for  the  2.0-MeV  radiation  simulation.  We  com¬ 
pared  simulation  results  to  four  of  the  characteristics  listed  in  the  key  of  Figure  4.2:  2DEG 
electron  mobility  /i 2deg>  threshold  voltage  Vth,  on  state  resistance  Ron,  and  transconduc¬ 
tance  Gm;Sat.  The  results  of  the  comparison  between  n 2deg  is  shown  in  Figure  4.1.  The 
comparison  between  the  model  outputs  and  experimental  results  are  laid  out  in  Figures  4.3, 
4.4,  and  4.5,  respectively. 
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Figure  4.3.  Comparison  Between  Model  Output  and  Experimental  Results 
for  Percent  Change  in  Vth  versus  Fluence 


Figure  4.4.  Comparison  Between  Model  Output  and  Experimental  Results 
for  Percent  Change  in  Ron  versus  Fluence 
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Figure  4.5.  Comparison  Between  Model  Output  and  Experimental  Results 
for  Percent  Change  in  Gmsat  versus  Fluence 


From  Figures  4.3,  4.4,  and  4.5,  as  stated  above,  the  model  overestimates  the  percent 
change  in  comparison  to  the  shifts  seen  in  the  physical  device.  This  could  be  because  of  the 
annealing  out  defects  in  the  physical  device.  The  physical  devices  were  not  tested  during 
the  radiation  event,  while  the  model  essentially  freezes  defects  in  place  and  does  not  allow 
annealing  out.  We  expect  the  trap  profile  to  begin  to  decrease  over  time  as  the  vacancies  and 
interstitials  anneal,  reducing  the  resultant  damage.  The  model  does  not  account  for  this. 

The  IV  characteristic  output  by  the  model  after  simulated  2.0-MeV  proton  irradiation 
to  a  fluence  of  6.0  x  1014  cm“2  are  shown  in  Figures  4.6,  4.7  and  4.8  along  with  the 
un-irradiated  model  outputs.  These  plots  clearly  show  the  threshold  voltage  shift  in  the 
device  towards  zero  at  both  100.0  mV  and  10.0-V  drain  voltage  values.  This  matches  the 
results  seen  in  physically  radiated  devices,  discussed  further  in  Section  4.3.  We  also  see  the 
current  saturate  at  a  much  lower  value  in  the  Id  versus  Vd  curve  shown  in  Figure  4.8,  which 
is  something  that  can  easily  affect  the  performance  of  a  component  utilizing  one  of  these 
devices. 
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Drain  Current  (A) 


Figure  4.6.  Vg  versus  Id  Comparison,  Vj  =  100.0  mV 


Figure  4.7.  Vg  versus  Id  Comparison,  Vd  =  10.0  V 
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Figure  4.8.  Id  versus  Vd  Comparison,  Vg  -  -3.0,  -2.0,  -1.0  and  0.0  V 


4.3  Trap  Occupancy  and  Modulation  of  the  Electric  Field 

In  order  to  determine  the  physical  origin  of  the  observed  shifts  in  the  device  characteristics, 
the  internal  structure  of  the  device  at  varying  fluence  was  examined.  A  review  of  the  struc¬ 
ture  files  produced  by  the  simulation  revealed  ionized  donor  and  acceptor  trap  occupancy 
mismatch  in  both  the  GaN  and  AlGaN  regions  of  the  device.  This  can  be  seen  in  Figures 
4.9  and  4.10  for  the  AlGaN  region  and  Figures  4.11  and  4.12  for  the  GaN  region  of  the 
device.  Both  acceptor  trap  defects  maintained  higher  ionization  levels  than  their  donor 
counterparts. 


39 


Figure  4.9.  Trap  and  Donor  Occupancy  Levels  (cm  3)  in  the  AlGaN  Barrier 
Layer  at  Fluence  level  of  10ucm-2 


Figure  4. 10.  Trap  and  Donor  Occupancy  Levels  (cm  3)  in  the  AlGaN  Barrier 
Layer  at  Fluence  level  of  6.0  x  1014cm-2 


40 


Microns 


Figure  4.11.  Trap  and  Donor  Occupancy  Levels  (cm  3)  in  the  GaN  Barrier 
Layer  at  Fluence  level  of  10ucm-2 


Microns 


Figure  4.12.  Trap  and  Donor  Occupancy  Levels  (cm  3)  in  the  GaN  Barrier 
Layer  at  Fluence  level  of  6.0  x  1014cm-2 
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The  trap  and  donor  ionization  difference  creates  a  build  up  of  negative  charge  in  the 
AlGaN  and  GaN  regions  of  the  device.  This  charge  buildup  screens  out  the  generation  of 
the  2DEG,  which  results  in  the  threshold  voltage  shift  towards  zero.  After  radiation,  less 
charge  is  needed  to  turn  the  device  off  as  the  2DEG  is  not  as  concentrated,  resulting  in 
decreased  current  for  a  given  bias  condition  and  decreased  transconductance. 

The  second  effect  is  seen  in  the  spread  of  the  electric  field,  as  shown  in  Figures  4.13 
and  4.14.  The  electric  field  under  the  gate  on  the  drain  side  of  the  device  spreads  out  over  a 
larger  area  at  higher  fluence  levels.  This  increase  in  width  decreases  the  field  and  increases 
the  breakdown  voltage  of  the  device.  The  effect  on  the  field  under  the  drain  edge  of  the  gate 
can  be  seen  more  quantitatively  in  Figure  4.15,  where  the  field  under  the  drain  side  of  the 
gate  is  shown. 


109.6  110  110.4 

Microns 


Figure  4.13.  Electric  Field  Distribution  at  a  Fluence  of  1011  cm  2 
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Figure  4.14.  Electric  Field  Distribution  at  a  Fluence  of  6.0  x  1014cm  2 


Microns 


Figure  4.15.  Electric  Field  Distribution  at  Fluences  of  1011  cm  2  and  6.0  x 
1014cm~2 
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4.4  Model  Evaluation  at  Additional  Energy  Levels 

Using  the  additional  radiation  profiles  derived  from  SRIM,  we  predicted  the  effects  that 
might  be  seen  in  an  actual  radiation  event  at  5.0,  10.0,  20.0  and  40.0  MeV  to  similar  fluence 
levels.  The  comparison  between  the  percent  change  in  key  device  characteristics  and  proton 
energy  levels  can  be  seen  in  Figures  4.16,  4.17,  4.18  and  4.19.  The  idea  that  the  higher  the 
energy  level  of  the  proton  the  more  fluence  is  required  to  inflict  a  similar  amount  of  damage 
to  the  device  was  expected  to  hold;  however,  the  increase  in  fluence  needed  to  shift  device 
characteristics  by  a  similar  amount  decreases  as  particle  energy  level  increases.  This  is 
unexpected,  as  decreasing  particle  interaction  should  make  it  harder  to  damage  the  device. 
As  the  particle  energy  increases,  we  expect  to  see  fewer  interactions  with  lattice  atoms, 
but  when  interactions  occur,  the  results  will  be  more  extreme  than  at  lower  energy  levels. 
Again,  this  is  an  issue  of  transferred  energy.  A  40.0-MeV  proton  has  a  lot  more  energy  to 
transfer  to  a  PKA  than  a  5.0- MeV  proton.  Where  a  5.0- MeV  proton  may  knock  three  atoms 
from  their  original  location,  a  40.0-MeV  proton  can  knock  out  50.  Also,  if  an  event  occurs 
with  a  higher  energy  level  proton,  the  proton  now  has  less  energy  and  chances  of  interacting 
with  other  particles  increase. 


Figure  4.16.  Percent  Change  in  pi2deg  versus  Fluence  for  Varying  Proton 
Energies 


44 


Figure  4. 17.  Percent  Change  in  Vth  versus  Fluence  for  Varying  Proton  Energies 


Figure  4.18.  Percent  Change  in  Ron  versus  Fluence  for  Varying  Proton  Ener¬ 
gies 
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Figure  4.19.  Percent  Change  in  gm,sat  versus  Fluence  for  Varying  Proton 
Energies 
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CHAPTER  5: 
Conclusions 


In  this  study,  we  have  successfully  created  a  TCAD  model  that  utilizes  Monte  Carlo  sim¬ 
ulated  proton  irradiation  damage  levels  for  several  different  proton  energy  levels  while 
incorporating  a  customized  calculation  for  the  electron  mobility  in  the  2DEG  to  simulate 
high-dose  proton  irradiation  in  a  GaN-on-Si  HEMT. 

Generally,  the  model  over-predicted  the  damage  that  was  observed  in  previous  NPS 
and  NRL  experimental  work.  We  did  not  over-predict  the  mobility  shift  in  radiated  devices 
because  the  mobility  was  used  to  tune  the  coefficients  in  the  customized  mobility  calculation. 
One  issue  with  using  this  metric  to  tune  the  model  is  that  the  devices  that  were  physically 
tested  were  not  immediately  characterized  after  irradiation.  This  allows  the  device  time 
to  anneal  out  defects  that  formed  during  the  irradiation.  This  can  result  in  a  shift  back  to 
less  damage  in  the  device.  The  model  effectively  freezes  the  defects  in  place  in  the  lattice 
without  any  function  to  allow  for  them  to  anneal;  however,  the  model  provides  a  good 
baseline  model  for  the  exploration  of  radiation  damage  in  GaN-on-Si  HEMTs. 

Continued  research  into  GaN-on-Si  devices  is  vital  to  the  Department  of  Defense  (DoD) 
as  more  and  more  hardware  systems  adopt  these  HEMTs.  They  strengthen  our  equipment, 
but  the  most  important  part  of  strength  is  knowing  weaknesses.  This  research  increases  our 
to  understanding  of  these  powerful  devices  that  are  now  being  utilized  in  active  systems. 

5.0.1  Future  Work 

The  Silvaco  software  does  not  have  specialized  formulas  which  depend  upon  2DEG  mobility 
such  as  those  proposed  in  research  by  Lisesivdin  [25].  For  future  work  on  this  project,  more 
physically  exact  terms  for  mobility  should  be  incorporated,  such  as  the  effects  of  roughness 
scattering.  For  the  HEMTs  that  were  used  in  this  study,  prior  work  by  Greenlee  [32] 
documented  the  visible  shift  in  atoms  at  the  boundary  layer  between  the  AlGaN  and  GaN 
layers,  which  has  an  effect  on  the  2DEG  in  that  region  through  roughness  scattering. 
Additional  research  should  also  examine  various  ways  to  model  the  electron  mobility  shift 
caused  by  PCF  without  physically  altering  the  model.  More  research  is  also  needed  in 
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this  simulation  to  model  the  diffusion  of  defects  throughout  the  lattice  in  order  to  identify 
mechanisms  for  void  sinks  causing  defects  like  the  voiding  seen  under  the  gate  of  the 
device  by  Wade  [10].  Lastly,  work  by  Pearton  identified  different  defect  types  created 
by  different  energy  levels  of  proton  radiation.  Additional  research  should  be  devoted  to 
identifying  these  trap  levels  and  incorporating  them  into  models  for  specific  energy  levels  of 
proton  irradiation.  This  model  can  serve  as  a  valuable  tool  for  predicting  the  effects  in  the 
experiments  that  should  follow,  40.0-MeV  proton  irradiation  of  NRL  GaN-on-Si  HEMTs 
to  fluence  levels  greater  than  1015cm-2.  This  tool  will  also  prove  valuable  for  upcoming 
neutron  irradiation  testing. 
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APPENDIX 


A.l  SRIM  Simulation  Values 

Table  A.l.  Interstitial  and  Vacancy  Fitting  Functions  in  the  GaN  Region  of 
the  Device 


Atom 

Slope  (( defects /cm) /(ion/ jam )) 

T-intercept 

Gcij 

395.816 

-83.213 

Gciy 

648.499 

-198.918 

Nj 

402.650 

-35.729 

Nv 

841.037 

-244.639 

5.0-MeV  Protons 


Table  A. 2.  Interstitial  and  Vacancy  Fitting  Functions  in  the  AlGaN  Region  of 
Device 


Atom 

Slope  {(de  f  ect s / cm) / (ion / ium)) 

T-intercept 

Alj 

918.307 

-388.069 

Aly 

102.064 

-19.034 

Gcij 

4299.410 

-1816.631 

Gay 

2817.0904 

-1117.765 

Ni 

10080.524 

-4296.504 

Ny 

2837.467 

-1042.749 

5.0-MeV  Protons 
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Table  A. 3.  Interstitial  and  Vacancy  Fitting  Functions  in  the  GaN  Region  of 
the  Device 


Atom 

Slope  (( de  f  ect  s  /  cm)  /  (ion  /  [mi )) 

K-intcrccpt 

Gcij 

315.347 

-91.402 

Gciy 

310.406 

-89.853 

Nj 

437.557 

-123.811 

Nv 

432.481 

-123.114 

10.0-MeV  Protons 


Table  A.4.  Interstitial  and  Vacancy  Fitting  Functions  in  the  AlGaN  Region  of 
Device 


Atom 

Slope  ((defects /cm) /(ion/ /uni)) 

K-intcrccpt 

Alj 

473.033 

-193.481 

Aly 

453.97 

-185.344 

Gcij 

2358.752 

-983.149 

Gciy 

2034.923 

-840.816 

N[ 

3345.905 

-1365.073 

Ny 

2752.402 

-1105.875 

10.0-MeV  Protons 


Table  A. 5.  Interstitial  and  Vacancy  Fitting  Functions  in  the  GaN  Region  of 
the  Device 


Atom 

Slope  (( de  f  ect s / cm) / (ion / nm )) 

K-intercept 

Gcij 

0 

31.988 

Gciy 

0 

31.968 

Ni 

0 

45.575 

Ny 

0 

45.085 

20.0-MeV  Protons 
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Table  A. 6.  Interstitial  and  Vacancy  Fitting  Functions  in  the  AlGaN  Region  of 
Device 


Atom 

Slope  ((de  f  ect  s  /  cm)  /  (ion  /  /um)) 

T-intercept 

Alj 

276.451 

-112.589 

Aly 

400.732 

-166.478 

Gcij 

1258.994 

-520.669 

Gciy 

1550.988 

-647.579 

Nj 

949.152 

-352.061 

Ny 

1454.471 

-572.491 

20.0-MeV  Protons 


Table  A.7.  Interstitial  and  Vacancy  Fitting  Functions  in  the  GaN  Region  of 
the  Device 


Atom 

Slope  ((de  f  ect s / cm) / (ion / /um)) 

T-intercept 

Gcij 

0 

17.561 

Gciy 

0 

17.238 

Ni 

0 

25.968 

Ny 

0 

25.489 

40.0-MeV  Protons 


Table  A. 8.  Interstitial  and  Vacancy  Fitting  Functions  in  the  AlGaN  Region  of 
Device 


Atom 

Slope  (( de  f  ect s  1  cm)  1  (ion  1  um )) 

T-intercept 

Alj 

-79.039 

40.018 

Aly 

-128.416 

61.496 

Gcij 

681.631 

-282.472 

Gciy 

686.527 

-284.880 

Nj 

790.314 

-310.940 

Ny 

844.579 

-335.516 

40.0-MeV  Protons 
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A.2  SILVACO  Simulation  Values 


Table  A. 9.  Albrect  Model  Values.  Adapted  from  [33] 


an 

0 

Bn 

5 

Cn 

0.4 

No 

1017 

To 

300.0 

T\ 

1065.0 

Values  A,  B  and  C  were  calculated  based  on  //„  of  1265  from  [34]  with  A  constantly 
zero. 


A.3  Trap  Levels  in  GaN 


Table  A.  10.  Defects  Introduced  by  Various  Types  of  Irradiation  in  GaN. 
Adapted  from  [20] . 


Material 

type,  dopant 

concentration 

(cm~  3) 

Radiation 

type 

Carrier 

removal 

rate, 

(cm~  1) 

Defect  levels 

(activation  en¬ 
ergy  (eV)  from 

conduction 

C  or  Valence 

V),  defect 

type  (donor  D, 
acceptor  A) 

Defect  pro¬ 
duction  rate 

(cm“l) 

Proposed 

identity 

n-GaN,  1016 

Protons, 

260 

0.13,  C,  ? 

30 

ER1 

2  MeV 

0.16,  C,  ? 

400 

ER2 

0.2,  C,  A 

600 

ER2 

n-GaN,  1016 

Protons, 

100 

0.2,  C,  ? 

ER3 

150  keV 

0.25,  C,  ? 

Complex  of  ER3 

0.32,  C,  ? 

Complex  of  ER3 

0.45,  C,  ? 

Complex  of  ER3 

0.6,  C,  A 

0.8,  C,D 

Gat 
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Table  A.ll.  Stated  Energy  Levels  for  Defects  in  GaN.  Adapted  from  [20]. 


Atom 

Vacancy  or  Interstitial 

Acceptor  or  Donor 

Deep  or  Shallow 

Energy  Level 

Ga 

Interstitial 

Donor 

Deep 

Ec  -  0.8  eV 

Ga 

Vacancy 

Acceptor 

Deep 

Ev  +  l.0eV 

N 

Interstitial 

Acceptor 

Deep 

Ec-  l.OeV 

N 

Vacancy 

Donor 

Shallow 

Ec  -  0.06  eV 

A.4  Code 

A.4.1  Silvaco  Deckbuild  Code 

go  atlas  simflags="-P  12" 

set  radenergy  =  2 

set  fluence  =  6el4 

set  talgan  =  17.5e-3 

set  tgan  =  0.8 

set  talgan2  =  0.517+0.258 

set  tain  =  0.427 

set  tsi  =  350 

set  tni  =0.03 

set  taug  =0.3 

set  tausd  =0.33 

set  tsin  =  0.133 

set  tair  =  0 . 366 

set  lsin  =  .133 

set  lgate  =  7 
set  Is  =  100 
set  lsg  =  3 
set  lgd  =  13 
set  Id  =  100 

set  Tem  =  21 
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set  pz  =  7.2el2 


mesh  width=100 

x.mesh  location=0  spacing=$ls/5 

x.mesh  location=$ls  spacing=$ls/100 

x.mesh  location=$ls+0 . 5*$lsg  spacing=$lsg/10 

x.mesh  location=$ls+$lsg-0 . 1  spacing=$lgate/18 

x.mesh  location=$ls+$lsg  spacing=$lgate/70 

x.mesh  location=$ls+$lsg+0 . 5*$lgate  spacing=$lgate/10 

x.mesh  location=$ls+$lsg+$lgate  spacing=$lgate/700 

x.mesh  location=$ls+$lsg+$lgate+0 . 1  spacing=0.02 

x.mesh  location=$ls+$lsg+$lgate+0 . 5*$lgd  spacing=$lgd/10 

x.mesh  location=$ls+$lsg+$lgate+$lgd  spacing=$ld/100 

x. mesh  location=$ls+$lsg+$lgate+$lgd+$ld  spacing=$ld/5 

y. mesh  location=-l*($tni+$taug+$tsin)  spacing=$tsin/2 
y.mesh  location=-l*($tausd)  spacing=$tausd/10 
y.mesh  location=-l*($tni)  spacing=$tni/10 

y.mesh  location=0  spacing=$talgan/10 

y.mesh  location=0 . 99*$talgan  spacing=$talgan/100 

y.mesh  location=0 . 999*$talgan  spacing=$talgan/1000 

y.mesh  location=$talgan  spacing=$talgan/100 

y.mesh  location=$talgan+0 . 002*$tgan  spacing=$tgan/250 

y.mesh  location=$talgan+0 . 02*$tgan  spacing=$tgan/200 

y.mesh  location=$talgan+0 . 2*$tgan  spacing=$tgan/20 

y.mesh  location=$talgan+$tgan  spacing=$tgan/6 

y.mesh  location=$talgan+$tgan+$talgan2  spacing=$talgan2/5 

y.mesh  location=$talgan+$tgan+$talgan2+$taln  spacing=l 

y.mesh  location=$talgan+$tgan+$talgan2+$taln+$tsi  spacing=$tsi/4 

region  number=l  material=AlGaN  x.mole=0.27  x.min=0  x.max=$ls+$lsg+$lgate+$lgd+$ld 
y.min=0  y .max=$talgan 

region  number=2  material=GaN  substrate  x.min=0  x.max=$ls+$lsg  y .min=$talgan 
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y .max=$talgan+$tgan 

region  number=3  material=AlGaN  x.mole=®.7  y .min=$talgan+$tgan 
y .max=$talgan+$tgan+$talgan2 

region  number=4  material=AlN  x.min=®  x .max=$ls+$lsg+$lgate+$lgd+$ld 
y . min=$talgan+$tgan+$talgan2  y . max=$talgan+$tgan+$talgan2+$taln 
region  number=5  material=Si  x.min=®  x.max=$ls+$lsg+$lgate+$lgd+$ld 
y .min=$talgan+$tgan+$talgan2+$taln  y .max=$talgan+$tgan+$talgan2+$taln+$tsi 

region  number=6  material=gold  x.min=®  x.max=$ls  y .min=-l*$tausd  y.max=® 
region  number=7  material=nitride  x.min=$ls  x .max=$ls+$lsg  y .min=-l*$tsin  y.max=® 
region  number=7  material=nitride  x .min=$ls+$lsg-Q . 1  x .max=$ls+$lsg 
y .min=- 1*  ($tni+$taug+$tsin)  y . max=- l*$tsin 

region  number=7  material=nitride  x .min=$ls+$lsg+$lgate  x.max=$ls+$lsg+$lgate+® . 1 
y .min=- l*($tni+$taug+$tsin)  y .max=-l*$tsin 

region  number=7  material=nitride  x .min=$ls+$lsg  x .max=$ls+$lsg+$lgate 
y .min=- 1* ($tni+$taug+$tsin)  y . max=- 1* ($tni+$taug) 
region  number=7  material=nitride  x .min=$ls+$lsg+$lgate 

x . max=$ls+$lsg+$lgate+$lgd  y .min=-l*$tsin  y.max=® 

region  number=8  material=nickel  x .min=$ls+Slsg  x.max=$ls+$lsg+$lgate 

y . min=-l*$tni  y.max=® 

region  number=9  material=gold  x .min=$ls+$lsg  x .max=$ls+$lsg+$lgate 
y .min=- 1* ($tni+$taug)  y . max=- l*$tni 

region  number=ll  material=gold  x .min=$ls+$lsg+$lgate+$lgd 

x . max=$ls+$lsg+$lgate+$lgd+$ld  y .min=- l*$tausd  y.max=® 

region  number=12  material=GaN  x .min=$ls+$lsg  x .max=$ls+$lsg+$lgate 

y . min=$talgan  y . max=$talgan+$tgan 

region  number=13  material=GaN  x .min=$ls+$lsg+$lgate 

x . max=$ls+$lsg+$lgate+$lgd+$ld  y .min=$talgan  y .max=$talgan+$tgan 
region  number=14  material=air  x.min=$ls  x.max=$ls+$lsg-Q . 1 

y . min=- 1* ($tni+$taug+$tsin)  y .max=-l*$tsin 

region  number=14  material=air  x.min^Q  x.max=$ls  y .min=-l,vC$tsin+$tausd) 
y . max= - 1 * ( $  t  au  s  d) 

region  number=15  material=air  x .min=$ls+$lsg+$lgate+® . 1 
x .max=$ls+$lsg+$lgate+$lgd  y .min=-l* C$tni+$taug+$tsin)  y .max=-l*$tsin 
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region  number=15  material=air  x .min=$ls+$lsg+$lgate+$lgd 
x .max=$ls+$lsg+$lgate+$lgd+$ld  y .min=- 1* ($tsin+$tausd)  y .max=-l*($tausd) 


electrode  name=gate  x .min=$ls+$lsg  x .max=$ls+$lsg+$lgate  y .min=-l*$tausd  y.max=® 
electrode  name=source  x.min=®  x.max=$ls  y .min=-l*$tausd  y.max=® 
electrode  name=drain  x .min=$ls+$lsg+$lgate+$lgd  x .max=$ls+$lsg+$lgate+$lgd+$ld 
y.min=-l*$tausd  y.max=® 
electrode  name=substrate  bottom 

doping  n.type  conc=lel8  uniform  x.min=®  x.max=$ls  y.min=®  y .max=$talgan 
doping  n.type  conc=lel8  uniform  x .min=$ls+Slsg+$lgate+$lgd 
x .max=$ls+$lsg+$lgate+$lgd+$ld  y.min=®  y .max=$talgan 

doping  f . doping="GaN_N_I . lib"  trap  acceptor  e.level=l.®  degen. fac  =  2 
sign=le-15  sigp=le-15 

doping  f . doping="GaN_N_V. lib"  trap  donor  e.level=3.34  degen. fac  =  2 
sign=le-15  sigp=le-15 

doping  f . doping="GaN_Ga_V. lib"  trap  acceptor  e.level=2.4  degen. fac  =  2 
sign=le-15  sigp=le-15 

doping  f . doping="GaN_Ga_I . lib"  trap  donor  e.level=2.6  degen. fac  =  2 
sign=le-15  sigp=le-15 

doping  f . doping="AlGaN_N_I . lib"  trap  acceptor  e.level=l.®  degen. fac  =  2 
sign=le-15  sigp=le-15 

doping  f . doping="AlGaN_N_V. lib"  trap  donor  e.level=3.34  degen. fac  =  2 
sign=le-15  sigp=le-15 

doping  f . doping="AlGaN_Ga_V. lib"  trap  acceptor  e.level=2.4  degen. fac  =  2 
sign=le-15  sigp=le-15 

doping  f . doping="AlGaN_Ga_I . lib"  trap  donor  e.level=2.6  degen. fac  =  2 
sign=le-15  sigp=le-15 

#doping  f . doping="AlGaN_Al_V. lib"  trap  acceptor?  e.level=2.4?  degen. fac  =  2 
sign=le-15  sigp=le-15  region=5 

#doping  f . doping="AlGaN_Al_I . lib"  trap  donor?  e.level=2.6?  degen. fac  =  2 
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sign=le-15  sigp=le-15  region=5 


#Values  taken  from  8-3,  8-4  from  pg  534  of  atlas  users  manual  using  GaN/2  and 
AIN/2  values  from  CH6  of  the  users  manual.  Exponent  term  was  not  halved, 
material  region=l  affinity=3.6  tcon. power  tc.npow  =  .317  tc. const  =  1.381 
material  region=2  affinity=3 . 932  tcon. power  tc.npow  =  .28  tc. const  =  .65 
material  region=3  affinity=2 . 78  tcon. power  tc. const  =  1.5535  tc.npow  =  .688 
material  region=4  affinity=2 . 03  tcon. power  tc.npow  =  1.64  tc. const  =  1.425 
material  region=5  affinity=4 . 05  tcon. power  tc. const  =  1.48 
material  region=6  tcon. power  tc. const  =  2.0 
material  region=8  tcon. power  tc. const  =  1.58 

material  region=12  affinity=3 . 932  tcon. power  tc.npow  =  .28  tc. const  =  .65 
material  region=13  affinity=3 . 932  tcon. power  tc.npow  =  .28  tc. const  =  .65 

contact  name=gate  work=4.6  surf.rec 
contact  name=source  work=3.6  resistance=10 
contact  name=drain  work=3.6  resistance=10 

models  print  k.p  fermi  srh 
models  lat.temp  heat. full 

interface  charge=$pz  s.s  y .min=$talgan  y .max=$talgan 

thermcontact  number=l  elec.number=4  temperature=273 . 15+$Tem 

mobility  region=2  f . tofimun="MobO . lib" 
mobility  region=12  f . tofimun="Mobl . lib" 
mobility  region=13  f . tofimun="Mob2 . lib" 

output  con. band  val.band  charge  int. charge  polar. charge  e. mobility 
method  block  newton  itlimit=25  maxtraps=20  nblockit=100  carriers=2 
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A.4.2  Tofimun  Mobility  Code 


GaN  Region  0 

#include  <stdio.h> 
#include  <stdlib.h> 
#include  <math.h> 
#include  <ctype.h> 
#include  <malloc.h> 
#include  <string.h> 
#include  <template.h> 

/* 


*  ATLAS  Parser  Function  Template 

*  ATLAS  Version  5.21. l.C 

*  c  1993  -  2015  SILVACO  Inc. 

*  All  rights  reserved. 


*/ 

/* 

*  General  field  dependent  mobility  model  for  electrons. 

*  TOtal  Field  Mobility  (parallel  and  perpendicular  field  components) 

*  Statement:  MATERIAL/MOBILITY 

*  Parameter:  F. TOFIMUN 

*  Arguments : 

*  Eperp  [in]  -  perpendicular  electric  field  (V/cm) 

*  Na  [in]  -  acceptor  concentration  (/cmA3) 

*  Nd  [in]  -  donor  concentration  (/cmA3) 

*  nconc  [in]  -  electron  concentration  (/cmA3) 

*  Eparl  [in]  -  parallel  electric  field  (V/cm) 

*  TL  [in]  -  lattice  temperature  (K) 

*  xcomp  [in]  -  x-species  fraction  (0-1) 

*  ycomp  [in]  -  y-species  fraction  (0-1) 
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*  *mun  [return]  -  hole  mobility  (cmA2/Vs) 

*  *dmundep  [return]  -  derivative  of  *mun  wrt  Eperp 

*  *dmundEparl  [return]  -  derivative  of  *mun  wrt  Eparl 

*  *dmundl  [return]  -  derivative  of  *mun  wrt  TL 

*  *dmundn  [return]  -  derivative  of  *mun  wrt  nconc 


V 


int  tofimun(double  Eperp, double  Na, double  Nd, double  nconc, double  Eparl, 
double  TL, double  xcomp, double  ycomp, double  *mun, double  *dmundep, 
double  *dmundepar , double  *dmundl , double  *dmundn) 

{ 

double  A=® . 25*2 . 61e-4 ,  B=4*2.9e-4,  C=4*17©e-4,  NQn=lel7,  T®n=30®.Q,  tln=l®65.0, 
alpha=Q.4,beta  =  1 . ® , Ecrit=2®e3  ,N1=2.®  ,Vsat=3e7; 
double  invmu , DD® , DD 1 ,DD2 , invmuALB ,h, g ,hPrime , gPrime ,MuField; 
double  N_imp=Na+Nd; 

invmuALB=A*(N_imp/N®n)*pow(TL/T®n, -1 . 5)+B*powC(TL/T®n) , 1 . 5)+C*pow((exp(tln/TL)-l) , -1) ; 
DD®=-1 . 5*A*(N_imp/N®n)*pow(TL/T®n, -2.5); 

DD 1=B* C 1 ■ 5)*pow ( (TL/TQn) , Q . 5) * ( 1/T®n) ; 

DD2=C*tln*pow((exp(tln/TL)-l) , -2)*(pow(TL, -2)) ; 

invmu= invmuALB ; 

h=l+pow((Eparl/Ecrit) ,N1) ; 

hPrime=(Nl/pow(Ecrit ,Nl))*pow (Eparl , CN1-1)) ; 
g=(l/invmu)+Vsat*(powCEparl , (Nl-l))/pow(Ecrit ,N1)) ; 
gPrime=(Vsat/powCEcrit  ,N1))  *CN1-1) -'powCEparl ,  (Nl-2))  ; 

MuField=g/h; 

,vmun  =  MuField; 

•'dmundepar  =  0 ; 

/ / (gPr ime/h) - ( 1/invmu) * (hPr ime/h) ; 

*dmundep=® . ® ; 
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* dmundl  = (DD0+DD 1+DD2) / (invmu*invmu) /h ; 
•'dmundn  =  0/(invmu*invmu)/h; 

return(0) ; 


GaN  Region  1 

Code  is  the  same  as  seen  in  GaN  Region  0  with  the  constant  values  changed.  Only  the 
constant  values  have  been  listed  here. 

double  A=0 . 25*2 . 61e-4 ,  B=1 . 35*2 .9e-4 ,  C=1 . 35*17Qe-4 ,  N0n=lel7,  T0n=300.0,  tln=1065.0, 
alpha=@.4,beta  =  1.0,  Ecrit=20e3,  Nl=2,  Vsat=3e7; 


GaN  Region  2 

Code  is  the  same  as  seen  in  GaN  Region  0  with  the  constant  values  changed.  Only  the 
constant  values  have  been  listed  here. 

double  A=@.25*2.61e-4,  B=2*2.9e-4,  C=2*170e-4,  N0n=lel7,  TOn=300.0,  tln=1065.0, 
alpha=@.4,beta  =  1 . 0 , Ecrit=20e3  ,N1=2.0  ,Vsat=3e7; 


A.4.3  Vacancy  /  Interstitial  Formation  Code 

AlGaN  A1  Interstitial  Code 

#include  <stdio.h> 

#include  <stdlib.h> 

#include  <math.h> 

#include  <ctype.h> 

#include  <malloc.h> 

#include  <string.h> 

#include  <template.h> 

/* 
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*  ATLAS  Parser  Function  Template 

*  ATLAS  Version  5.21. l.C 

*  c  1993  -  2015  SILVACO  Inc. 

*  All  rights  reserved. 


V 

/* 

*  Position  dependent  net  doping. 

*  Statement :  DOPING 

*  Parameter:  F. DOPING 

*  Arguments : 

*  x  location  x  (microns) 

*  y  location  y  (microns) 

*  *nnet  net  doping  concentration  (per  cc) 

*/ 

int  doping (double  x, double  y, double  *nnet) 

{ 

double  talgan  =  0.0175; 
double  tgan  =  0.8; 
if (y<0) { 

*nnet  =  0 ; 

} 

else  if(y>=@  &&  y<=talgan){ 
double  flu=8e+14; 

*nnet=  flu*5®15 . 561905*(y+0 . 43) -2143 . 1099108 ; 

} 

else{ 

,vnnet=0 ; 

} 

return(0)  ; 

} 
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